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INTRODUCTION 


This report provides the experimental results obtained by the use 
of ultrasonic techniques as a non destructive evaluation tool for fiber 
reinforced composite materials. The report is divided into two parts, 
because the research was done in two separate periods, in two indepen- 
dent locations, and using two different pulsing systems. 

The first part covers the ultrasonic evaluation of a "0" degree 
graphite-epoxy fiber reinforced composite under tensile loading. The 
research proved that wave speeds are frequency dependent, that attenuation 
is strongly dependent on tensile stress, and that wave speed is weakly 
dependent on tensile stress. 

The second part is devoted to understanding the dispersion mechanism 
in thin graphite-epoxy fiber reinforced composite. This dispersion mecha- 
nism has to be fully examined and understood in this material in order to 
sort out the relationships existing among stress, strain, damage, and 
ultrasonic signals. 
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ABSTRACT 


An Acoustical -Ultrasonic Technique was used to demonstrate that 
relationships exist between changes in attenuation of stress waves and 
tensile stress for an 8-ply "0" degree graphite-epoxy fiber reinforced 
composite. All tests were conducted in the linear range of the material 
for which no mechanical or macroscopic damage was evident. Changes in 
attenuation were measured as a function of tensile stress in the frequency 
domain and in the time domain. Stress wave propagation in these specimens 
was dispersive, i.e., the wave speed depends on frequency. Wave speeds 
varied from 267400 cm/sec to 680000 cm/sec as the frequency of the signal 
was varied from 150 KHZ to 1.9 MHZ which strongly suggests that flexural/ 
Lamb wave modes of propagation exist. The magnitude of the attenuation 
changes depended strongly on tensile stress. It was further observed that 
the wave speeds increased slightly for all tested frequencies as the stress 
was increased. 
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CHAPTER I 
INTRODUCTION 


The first practical application of ultrasonics for test metal bodies 

for internal defects was done by the Russian Sokolov^ in 1934. He was 

using a through transmission technique, transmitting the sound into the 

body under test, utilizing a piezo-electrical crystal, and was detecting 

this sound on the opposite side of the body, adopting a second similar 

piezo-electrical crystal. Independently in .this country as well as in 

England ultrasonic development started about 1940. 

After World War II the measurement of stress wave attenuation by means 

of ultrasonic nondestructive evaluation (NDE) techniques evolved gradually 

to a more complete and desirable state. In the 50 ‘s NDT techniques were 

used for flaw detection and for measuring the thickness of material where 

only one side was accessible. Also automatic monitors were adopted in the 
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instrumentation of ultrasonic testing. In 1954, Firestone developed the 

3 

principle for the impulse echo devices. In 1956, Truell and Hikata esta- 
blished the relationship between attenuation in aluminum and the number of 
fatigue cycles. In 1960 ultrasonic determination of elastic constants of 

4 

metals at elevated temperature was achieved by Levitt and Martin . Ultra- 
sonic pulse-echo techniques have been used by evaluating thickness, bonding, 
and corrosion . 

Recently, (NDE) techniques have gone further than detecting flaws; Vary 
et.al.^"^^ have demonstrated that a correlation does exist between ultra- 
sonic attenuation and mechanical and strength properties of metals, ceramics, 

IP 1 

and fiber composites. Also Williams et.al. ’ have shown that ultrasonic 
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attenuation is an indicator of a fatigue life of graphite-epoxy composites. 
Ultrasonic measurements in the area of composite material, up to now, have 
provided a nondestructive means for measuring variations in strength related 
properties, predicting interlaminar shear strength of fiber composite 
laminates, and ranking composite structures according to strength. 

This research will investigate how well (NDE) techniques describe stress 
related effects. An ultrasonic attenuation study, of a stress wave propa- 
gating through a "0" degree graphite-epoxy fiber reinforced composite specimen 
subjected to tensile loads, was performed. This paper will provide not only 
qualitative data, but also quantitative data, showing the relationships that 
exist between changes in attenuation of stress waves and tensile stress for 
the 8-ply "0" degree specimen under consideration. Tensile testing was done 
in the linear range of the material, where linearity is defined by constant 
slope of the stress-strain curve in the range of the test conducted. Wave 
speeds were also measured showing that stress wave propagation is dispersive. 

This thesis will demonstrate that: 

1. Dispersion does occur, i.e., wave speeds are frequency 
dependent. 

2. The attenuation is strongly dependent on tensile stress. 

The wave speed is weakly dependent on tensile stress. 


3. 


3 


CHAPTER II 

EXPERIMENTAL PROCEDURE 


The specimen used was an eight-ply "0" degree graphite-epoxy fiber 
reinforced composite. An AS graphite fiber and a PR-288 exposy resin 
have been used to construct an AS/PR-288 preimpregnated fiber-resin ply. 
The fibers were coated with a polyvinyl alcohol to enhance fiber-matrix 
interface strength properties. A schematic of the unidirectional compo- 
site is shown in Figure 1. Physical properties of the specimen are given 
in Table I below:® 

TABLE I - DESCRIPTION OF THE TENSILE SPECIMEN LVll 


Ply angle, deg® 

[0]8 

Ultimate strength, GP, 

a 

1.24 

Fiber coating 

PVA 

Tensile modulus E, GP 

a 

no 

Gage length, cm 

12.7 

Poisson's ratio. 

.283 

Thickness t, cm 

.1384 

Fiber fraction 

69.2 

3 

Density p, g/cm 

1.584 

Stress wave factor NSW 

.853 


® Fiber angles relative to tensile axis in AS graphite, 

PR-288 epoxy composite 

f^ determine by HNOj digestion (ASTM-D3171) 

2.1 APPARATUS: 

A schematic of specimen and transducers in ultrasonic through trans- 
mission testing was shown in Figure 2. 



(Through thickness 
direction) 
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► 

(Transverse 

Direction) 


(0*; Longitudinal direction) 


FIG. 1 — PRINCIPAL DIRECnCWS OF 
"0" DEGREE CSAPHITE 
EPOXY SPECIMEN 
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Panametric transducers as well as Dunegan transducers were used in 
this experiment. Their properties are tabulated in Table II. 


TABLE II - THE PROPERTIES OF THE TRANSDUCERS 


Frequency, 

Model 


MHZ 2.25 

Panametri c 
V314 


1.0 

Panametric 


A306 (-) 


.438 

Dunegan 

S750B 


.226 

Dunegan 

S140B 


Transmits 
& Receives 


Longitudinal 

Waves 


Longitudinal 

Waves 


Longitudinal 

Waves 


Longitudinal 

Waves 


Type 


Broad-band Broad-band 


Resonant 


Resonant 


Aluminum semi cylindrical buffers of 3.157 mm (.125 in) diameter were 
attached to the transducers in order to: 

a. Provide a line contact of length 12.7 mm (.5 in) which 
the width of the specimen. 

b. Prevent bending from occurring as pressure was applied 
to the transducers by the transducer holders. 

c. Place the sending and the receiving transducer at specified 
location. 


Strain gauges were mounted on the specimen on both faces to measure 

the strain when tensile testing was performed. Strain gauges were 0° and 

90® orthogonal back-to-back. Also strain gauges were installed onto the 

transducer holders to measure the pressure applied to the transducers. 

13 

Williams, et. al. defined and measured a "saturation pressure" which is 
required for reproducible coupling of transducer to specimen; this pressure 
was monitored by the strain-gauged transducer holders. 

An ultragel II ultrasonic couplant was used between the aluminum buffer 
and the specimen. A very thin coating of couplant was used where effects of 
thickness of the couplant were neglected. A felt material was used on the 
backside of the transducer holder to insure that energy did not enter the 
transducer holder from the specimen. 
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2.2 EXPERIMENTAL SYSTEM 

The experimental system is shown in Figure 3 and the description of 
the main equipment is tabulated in Table III. 


TABLE III - DESCRIPTION OF THE MAIN EQUIPMENT 


IDENTIFICATION 

Tensile 

Parts: a - Honeywell machine 

b - Riehle machine 

Ultrasonic 

Parts: a - Spectrum analyser 

7L12 

b - Panametric Pulser 
5052 P/R 

c - Krohn-hite filter 


d - Dunegan totalizer 
301 

e - Oscilloscope 
7B92 


FUNCTION 

a - Strain and load indicator 
b - Tensile machine 


a “ Outputs the frequency spectrum 
of the transmitted pulse 

b - Provides a voltage pulse to 
sending transducer 

c - Provides band-pass filtering of 
transmitted pulse 

d - Amplifies and band-pass filters 
transmitted pulse 

3 - Records the input pulse and the 
transmitted pulse in real time 


2.3 PROCEDURE 

An ultrasonic wave was transmitted into the specimen by a sending 
transducer. Two inches away in the longitudinal direction, and on the 
same face, a receiving transducer was placed to sense the stress wave 
energy traversing the specimen. The receiving transducer output was 
displayed on the oscilloscope for real time analysis and on the spectrum 
analyser for frequency analysis. 

The specimen, mounted in a tensile machine where care was taken to 
avoid off-axis loading or twisting, was pulled at an initial loading rate 
of approximately 17pm/sec (.04 in/min) to set the grips. After this 




SCOPE 



SYSTEM FOR VELOCITY AND 
ATTENUATION MEASUREMENTS 


FIG. 3 — EXPERIMENTAL SYSTEM -1 
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initial loading to about 445N (100 Ibf) which was the reference load 
for attenuation measurements and change in attenuation measurements, 
further loading was done at mean rate of 4.2 jjm/sec (.01 in/min). 

The load was increased to 3560N (800 Ibf), 6675N (1500 Ibf), and in 
some experiments to 13350N (3000 Ibf) successively. Then the load was 
decreased to 6675N, 3560N, and 445N. At the same time strain-load data 
was generated, photographs of the stress wave output and the frequency 
output were taken. 

The loading-unloading cycle was repeated five to seven times for 
each setting of independent variables (Table IV), which usually occurred 
within a two hour period. The number of cycles that were required was 
dictated by the reproducability of the attenuation changes. 

Typical output signals are shown in Figure 4 and Figure 5 for an 
entire loading-unloading cycle. 


TABLE IV - SETTING OF INDEPENDENT VARIABLES 


Frequency® 

Sending^ 

Transducer 

n . . b 

Receiving 

Transducer 

Number® 
of Tests 

Filter‘d 

Setting 

1.9® 

2.250 

2.250 

2 

2 MHZ/2MHZ 

.850 

1.000 

.439 

2 

1 MHZ/1 MHZ 

.250 

2.250 

1.000 

2 

.3MHZ/.4MHZ 

.180 

1.000 

.439 

2 

.2MHZ/.2MHZ 

.150 

.439 

.226 

2 

.IMHZ/.IMHZ 


a - From spectro-analysis of the output signal (in MHZ) 

b - Frequency of transducer type (in MHZ) 

c - Number of independent tests performed in different days 

d - Krohn-hite filter was used passing everything > to x and 
•< to y in the setting x/y. 

e - Also 1.7 MHZ was displayed in the same output with 1.9 MHZ. 
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FIG. 4 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of the 150 KHZ in 
the Frequency Domain 
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n 



d. At 361 7N 



e. At 525N 


FIG. 4 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of 150 KHZ in 

the Frequency Domain 
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a. At 507N 


b. At 3658N 



FIG. 5 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of the 150 KHZ 

in the Time Domaine 



flOOnV 


StOnS 



FIG. 5 - ULTRASONIC PULSE ATTENUATION VERSUS 

the Time Domain 
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d. At 361 7N 


e. At 525N 


OF 150 KHZ in 
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CHAPTER III 
RESULTS 

In this experimental research many observations were made concerning 
how various parameters effected the ultrasonic output; only a few were 
examined and are presented in this paper because of time, equipment, and 
experimental system limitations. This chapter will present the disper- 
sion phenomenon and the wave type, the attenuation as related to tensile 
load, and the wave speed as related to tensile load; the specimen was a 
"0" degree graphite-epoxy fiber reinforced composite. Also reproducibility 
of results will be treated at the end of this chapter. 

3.1 DISPERSION AND WAVE TYPE 

Dispersion is understood to mean that wave speeds are frequency 

dependent. From the data collected and the measurements taken, a graph 

was plotted showing speed versus frequency for the specimen under study 

in Figure 6. The wave speed was determined by measuring the delay time 

between the output wave front and input wave front over fixed distance. 

The graph shows that the wave velocity along the longitudinal fiber is 

frequency dependent. Wave speeds varied from 267400 cm/sec to 680000 cm/sec 

as the frequency of the signal was varied from 150 KHZ to 1.9 MHZ. 

There are actually two types of dispersions. First is viscoelastic 

dispersion due to viscoelastic properties of the material. Second is a 

13 

geometric dispersion due to the dimensions of the material. Williams 
demonstrated that there is no viscoelastic dispersion in the graphite-epoxy 
composite for longitudinal and shear wave velocities in the range tested: 
Shear wave velocities tested in the range of 0.48 to 3.0 MHZ and longitudinal 



850000 

680 000 - 

628000 - • 
502 500 - • 

357200 - • 

267 400 - • 


180 

Lli. I I 

150 250 850 

Frequency f, KHZ -*• 


FIG. 6 WAVE SPEED VERSUS FREQUENCY 
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wave velocities tested in the range of 0.225 MHZ to 5 MHZ. The frequency 

range of results presented in this paper is 0.15 to 1.9 MHZ, which is 

very close to the range of Williams' result. 

For the wave propagating in the x-j direction and having a particle 

13 

motion in the direction, William measured the longitudinal wave 

c 

speed and found it to be 8.5 to 10 cm/sec. This wave speed is 
plotted in Figure 6 and is used as a comparison since the material used 
in William's research is very similar to the material used in the 
research reported here. 

By comparing Cj^ to the velocities obtained in this experiment as 

shown in Figure 6, one can see that the wave velocities measured were 

not longitudinal wave velocities. The sending transducer inputs a 

longitudinal wave perpendicular to specimen which, because of many 

reflections, becomes a wave propagating along the length of the specimen. 

During this propagation shear waves are formed due to the reflections 

and interact with longitudinal waves. This interaction of waves propa- 

14 

gating along a thin specimen produces Lamb waves . Dispersion does 

15 

occur in the propagation of flexural /Lamb waves . 

Lamb wave propagation may exist in a multimode state at a given frequency 
producing a finite number of wave speeds to exist simultaneously. Because of 
experimental limitations, this research could not define which mode (Figure 6) 
and how many modes were superimposed upon each other, but this paper recognized 


II I 


I I III 


I I HI 


II 
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the possible existence of different modes simultaneously in the stress 
wave output. 

3.2 ATTENUATION IS STRONGLY DEPENDENT ON TENSILE STRESS 

The relationship between attenuation and tensile stress was investi- 
gated by measuring the pulse amplitude change in the time domain as well 
as in the frequency domain. The tensile load applied to the specimen was 
in the linear range of the material where linearity was based on the constant 
slope of the stress-strain curve shown in Figure 7. 

In order to understand the analysis adopted in this chapter, a graph 
is shown in Figure 9 which refers to general but realistic pulse. Ar.d this 
graph is complemented by the following explanation: 

The percent change in amplitude in the time domain (Figure 8) was cal- 
culated in the following manner — the peak to peak amplitude A^ for a 
spike was measured at every specified load P^; then the difference between 
Ag and the reference peak to peak amplitude Aj^ at the reference load Pj^ 

(first 445N load reached in the cycling process) was multiplied by the 
inverse of Aj^ and by one hundred which gave the percent change in amplitude 
at Pg. Most of the time though, this analysis was applied to n number of 
spikes in a stress wave output at Pg: 


% amplitude change at Pg 



n 

l 

i=l 



X 100 


Equation 1 
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STRESS WAVE OUTPUT 
IN THE TIME DOMAIN 



FREQUENCY OUTPUT 
IN THE FREQUENCY 
DOMAIN 


FIG. 8 SCHEMATIC OF Ag AND A^ IN THE 

TIME DOMAIN AND IN THE FREQUENCY 
DOMAIN 
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where Pg = Specified load 445N, 3560N, or 6675N... 

i =1, — , n. Where n is the number of spikes. 

= Vertical measurement from peak to peak of a specified 
^ spike when a specified load is applied. 

An = Vertical measurement from peak to peak of a reference 
spike when the reference load is applied. 

The number of spikes chosen was dictated by observing which group of 
spikes behaved in a similar way and were of the same frequency. 

The same analysis procedure was applied to the frequency domain 
response (Figure 8), the only difference was that only one spike existed 
at each frequency, and A^ or A^ would be the amplitude. Consequently, 
Equation One becomes: 

% amplitude change at = [(A^ - x 100 Equation 2 

Two different phenomenon were observed as the load was cycled. First, 
during the first few cycles of loading the relationships between amplitude 
changes and load shifted upward from line A of Figure 9. After 3 or 4 cycles 
the shifting ceased and the data followed line Cj the data being contained in 
the scatter bandwidth B. Second, the percent amplitude change was behaving 
linearly, increasing or decreasing depending on the frequency of the output 
signal, when tensile loads were increasing. Figure 9. 
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FIG. 9 ~ GSy^H DESCRIBING ANALYSIS PROCEDURE 

T = Total change A * First cycle line 

B « Band width C = Stable cycle line 

S = Shift 
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3.2.1 FREQUENCY DOMAIN RESULTS 

Figure 10 shows the relationship between the percent change in 
amplitude versus tensile load for the six frequencies under considera- 
tion in the frequency domain. It is approximately a linear relationship, 
increasing with load for the 150, 180, 250 and 1900 KHZ, and decreasing 
when load is increasing for the 850 and 1700 KHZ, which demonstrates that 
attentuation is not only dependent on tensile load, but it is also frequency 
dependent. Table V describes shift, B/2 and T for every frequency under 
study. 

TABLE V - SHIFT, BANDWIDTH, AND TOTAL CHANGE FOR EVERY FREQUENCY 

UNDER CONSIDERATION 

^KHZ S B/2 T 


150 

13 

3 

12 

180 

16 

5 

8 

250 

13 

4 

20 

850 

28 

7 

28 

1700 

48 

5 

20 

1900 

20 

8 

25 


3.2.2 TIME DOMAIN RESULTS 

Figure 11 shows the relationship between the percent change in amplitude 
versus tensile load for the frequencies under consideration. Table VI shows 
B,S, and T for the time domain results which may be compared with the frequency 
domain results of Table V. 


TABLE VI 

- BEHAVIOR OF FREQUENCIES 

IN THE TIME DOMAIN 


^KHZ 

S 

B/2 

T 

150*^ 

14 

3 

14 

150a 

0 

6 

10 

180 

13 

4 

5 

250 

11 

5 

18 

850 

14 

6 

16 

1900,1700 

20 

3 

8 
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FIG. 10 — PERCENT CHANGE IN AMPLITUDE 
VERSUS TENSILE LOAD IN THE 
FREQUENCY DOMAIN 

- Note all lines in this graph are of "A" type (See Fig. 9) 
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PERCENT AMPLITUDE CHANGE 


24 



FIG. 11 ~ PERCENT CHANGE IN AMPLITUDE 
VERSUS TENSILE LOAD IN THE 
TIME DOMAIN 

- Note all lines in this graph are of "A" type (See Fig. 9) 
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The only exception was the 150 KHZ case where in experiment a, a definite 
amplitude increase with load was observed, and in experiment b, a definite 
decrease in amplitude was seen when load was increased. Photographs for 
the same cycle in both domains are given in Figure 12 and in Figure 13 for 
the second experiment which could be compared to the photographs in Figure 
4 and Figure 5 in Chapter 2. 

3.2.3 COMPARISON OF RESULTS BETWEEN FREQUENCY DOMAIN AND TIME DOMAIN 

Comparing results from both domains has helped not only to show that 
both domains behaved in the same way in terms of total energy analysis, but 
also has helped to observe the presence of superposition of different 
frequencies as well as of different modes in the stress wave output in 
the time domain. 

3.2. 3.1 SUPERPOSITION OF DIFFERENT FREQUENCIES 

An extensive calculation on individual spikes was done in the case of 
the 1.9 and 1.7 MHZ experiment, and showed that some spikes were increasing 
with load while some others were decreasing. This behavior was analysed in 
comparison with the behavior of the two frequencies 1.9 and 1.7 MHZ in the 
frequency domain, and results were tabulated in Table VII and Table VIII. 
Also photographs were presented for this case in Figure 14 and Figure 15. 
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b. At 3587^ 


c. At 6786N 




FIG. 12 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD Of The 150 KHZ 

in the Frequency Domain 
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FIG. 12 - ULTRASONIC PULSE ATTENUATION V! 

in the Frequency 
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b. At 3587N 




FIG. 13 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of the 150 KHZ 

in the Time Domain 
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At 3565N 
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FIG. 13 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD OF 150 KHZ in 

the Time Domain 
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TABLE VII 

- BEHAVIOR 

OF THE 1 . 

7 MHZ SPIKES IN 

THE TIME DOMAIN 




THE 

SPIKE NUMBER 


L 


21st 

22nd 

23rd 

24th 

0 


% 

% 

% 

% 

A 






D 

445 

0.00 

+10.71 

+35.29 

+20.00 

A 






P 

3560 

-17.24 

-10.71 

+ 5.88 

-13.33 

P 






L 

I 

E 

6675 

-27.59 

-28.57 

-11.76 

-40.00 

3560 

-13.59 

-10.71 

+17.65 

-13.33 

D 

1 

445 

+ 6.90 

+10.71 

+64.71 

+53.33 


N 


TABLE VIII-BEHAVIOR OF THE 1.9 MHZ SPIKES IN THE TIME DOf^IN 


L 

0 

A 

D 

A 

P 

P 

L 

I 

E 

D 


445 

3560 

6675 

3560. 


THE SPIKE NUMBER 


8th 

9th 

10th 

% 

% 

% 

-25.64 

17.39 

11.11 

-12.82 

26.09 

27.78 

+ 2.56 

43.48 

27.94 

-10.26 

30.43 

22.06 

-30.77 

17.39 

11.76 


11th 12th 13th 

% % % 

5.88 2.60 3.64 

17.65 14.29 21.82 

27.94 19.48 30.91 

15.58 15.58 25.45 

7.79 12.73 


14t 

% 

4.6 

23.2 

41.8 

27.9 

16.2 


N 


445 


7.79 
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a. 



b. 



c. 


At 445N 


At 3560N 


At 6675N 


FIG. 14 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of the 1.7 and 1.9 MHZ in the 

Frequency Domain 
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d. At 3560^ 



FIG. 14 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of the 1.7 and 1.9 MHZ 

in the Frequency Domain 
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a. At 445N 



b. At 3560N 
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FIG. 15 - ULTRASONIC PULSE ATTENUATION VERSUS 

MHZ in the Time Domain 


d. At 3560N 


e. At 445N 


of the 1 .7 and 1 .9 
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3. 2. 3. 2 SUPERPOSITION OF DIFFERENT MODES 

Even when spectro-analysis results indicated a very narrow band of 
frequencies in the output signal, the time analysis on an individual spike 
base indicated the presence of different modes: Some spikes were behaving 

exactly as the single spike in the frequency domain, some others were behaving 
differently and in the rest of the spikes no definite behavior was observed. 
And as mentioned previously in Section 2.1, lamb waves may exist in multimode 
states at a given frequency. For better understanding of what was said pre- 
viously in this section one could follow the behavior of some spikes in 
Figure 5 and compare it to the result in Figure 10 or Figure 11. 

3.3 WAVE SPEEDS ARE WEAKLY DEPENDENT ON TENSILE LOAD 

From the data collected, the time shift for the first early spikes was 
measured at each load having 445N as our reference load. Thereafter a per- 
cent of change in wave speed was calculated and is presented in Table X below. 

TABLE IX - PERCENT OF CHANGE IN WAVE SPEED 


F 

445 

L0AD,N 

3560 

6675 

R 

q 150000 

0 

1.87 

3.20 

g 180000 

0 

1.03 

1.77 

^ 250000 

0 

.72 

.84 

y 850000 

0 

1.02 

1.43 

5 1900000 

0 

1.70 

3.40 

When tensile loads were applied the distance between the two transducers 

increased and data showed 

that the 

corresponding delay time got 

smaller in 
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Table IX. The phenomenon described cannot be explained at this time using 
linear theory. 

3.4 REPRODUCIBILITY 

Reproducibility means repeatability of the changes that were observed 
and their stability when cycling was occurring. As was mentioned earlier 
only two formal experiments for each frequency were recorded. In all experi- 
ments except for the 150 KHZ experiment, repeatability was achieved after 
the cycling process had stabilized. As an example the 180 KHZ first experi- 
ment is presented in Figure 16 and Figure 17, and the second experiment in 
Figure 18 and Figure 19, where one could observe and follow the changes in 
both experiments in both domains and perceive how well they agreed. 

Reproduciblity of the signal itself was not a primary concern because 
of the experimental system limitations: 

a - Coupling thickness was controlled only by the line 
pressure of the aluminum buffers 

b - Riehle tensile machine could not be controlled with 
enough precision. 

c - Distance between transducers as well as placing those 
transducers on the specimen did cause differences in 
the output signal itself. 

Even though limitations existed, good case was taken to reproduce almost 
the same signal as presented in Figure 17 and Figure 19. 
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a. At 458N 


b. At 3630N 


c. At 6786N 




FIG. 16 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of the 180 KHZ in 

the Frequency Domain 
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FIG. 16 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of 180 KHZ in the 

Frequency Domain 



b. At 3630N 


c. At 6786N 
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FIG. 17 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of 180 KHZ in the 

Time Domain 


1 





FIG. 18 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of the 180 KHZ in 

the Frequency Domain 
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FIG. 18 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of 180 KHZ in the 

Frequency Domain 







c. At 6888N 


FIG. 19 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of the 180 KHZ in 

the Time Domain 
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FIG. 19 - ULTRASONIC PULSE ATTENUATION VERSUS LOAD of 180 KHZ in the 

Time Domain 
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CHAPTER IV 
CONCLUSION 


It has been shown that dispersion did occur, i.e., wave speeds 
were frequency dependent. Moreover, it was found that flexural /Lamb 
waves existed, and a superposition of different modes existed in the 
stress wave output of the signal. 

It has been proven that attenuation depended strongly on tensile 
stress. Also it was observed that attenuation was frequency dependent. 

Finally, it was demonstrated that wave speeds were weakly dependent 
on tensile load. 

Further research is needed using a pulsing system producing narrow 
band frequencies since all effects measured were frequency dependent. 



PART II 


LAMB WAVE PROPAGATION 
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CHAPTER 1 
INTRODUCTION 


The research conducted here is devoted to gaining an understanding 
of the geometric dispersion mechanisms in fiber reinforced composite 
specimens. Geometric dispersion means that wave speeds are frequency 
dependent because of geometric constraints and not because of material 
properties. A short review of wave propagation in solids will be presented 
as an introduction to the research results. 

The theory of elastic wave propagation in solids shows that two types 
of elastic waves may be propagated in an unbounded or extended isotropic 
medium - longitudinal waves where the particle motion is along the direc- 
tion of propagation, and transverse waves where the particle motion is 
normal to the direction of propagation or parallel to the wave front. Each 
of these two wave types has its characteristic wave speed and is non dispersive. 

In extended or unbounded anisotropic elastic materials three waves of 

different speeds may exist for wave propagation in an arbitrary direction. 

In general, the particle displacement has components along and normal to the 

direction of wave propagation for all three waves. These waves cannot be 

classified as longitudinal or transverse waves except in certain specified 

directions corresponding to axes of symmetry. The wave propagation is non 

1 3 

dispersive in linear elastic anisotropic materials. Williams et.al. 
has shown that dispersion does not occur in the graphite-epoxy materials 
he tested. 
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Elastic surface waves may exist where there is a bounding surface. 

This type of wave is created because the free boundary conditions must 
be satisfied. These surface waves travel without change in form; the 
particle motion is both vertical and horizontal with an elliptic tra- 
jectory. Moreover they are frequency independent or nondispersive. 

Lord Rayleigh^® has shown that their effect decreases rapidly with depth 
and that their velocity of propagation is smaller than that of transverse 
waves . 

In an isotropic plate of small thickness the analysis of wave propa- 
gation becomes more complicated. Wave propagation in plates or waveguides 
produces what are called Lamb (or plate) waves which exhibit geometric 
dispersion. Once the waves start bouncing back and forth between the 
bounding surfaces it is very difficult to trace the actual reflections. 

The general direction of propagation is, however, parallel to the bounding 
surfaces. The analysis of harmonic wave propagation in wave guides leads 
to a dispersion mechanism with modes of wave propagation where wavelength 
to thickness ratio, phase velocity, and group velocity are important quan- 
tities. 

When an incident longitudinal or transverse wave strikes a free surface 
two waves are generated on reflection, one is longitudinal, the other is 
shear wave. The interaction between longitudinal and shear waves in a 
relatively thin plate produces Lamb waves travelling parallel to the free 
surface. These waves in which the particle motion occurs both parallel to 
the wave propagation and perpendicular to the freesurface, are of two families 
of modes, symmetric and antisymmetric (see Figure 20). The symmetric modes 
have particle motions which are symmetric with respect to the midplane, i.e. 
for (x,y,z) and (x,y,-z) the longitudinal displacements are equal, but the 
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a. antisymmetrical 
lamb wave 



b. symmetrical 
lamb wave 


Fig. 20 


Symmetrical and antisymmetrical lamb waves. 



particle motion 
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transverse displacements are equal in magnitude and opposite in direction. 

The antisymmetric modes have particle motions which are antisymmetric with 
respect to the midplane, i.e. for (x,y,z) and (x,y,-z) the transverse dis- 
placements are equal in magnitude and opposite in direction. Infinite 
number of modes of propagation exist in a plate according to Horace Lamb^^, 

TO 

this was confirmed experimentally by D. C. Worlton for aluminum and 
zirconium materials. 

Rayleigh waves or Lamb waves may exist on the surface of a wave guide 

depending on the thickness of the wave guide and the frequency of excitation. 

As the product of frequency and thickness gets large, the boundary conditions, 

in effect, decouple and the Lamb waves become Rayleigh waves as shown in 

19 

Figure 21. Meitzler has demonstrated experimentally that different modes 

of Lamb wave propagation exhibit different attenuation behavior because of 

the strong frequency and thickness dependencies. 

20 

Habeger, Mann, and Baum have studied Lamb wave propagation in paper, 
an orthotropic material. They have developed the dispersion curves for the 
symmetric and the antisymmetric cases. Their analysis may be applied directly 
to unidirectional fiber reinforced composites. 

The review of literature on wave propagation in solids and the results 
presented in Part I of this report clearly demonstrate that wave propagation 
in thin plates is in the form of Lamb waves. The research reported in this 
part will present experimental evidence of Lamb wave propagation and related 
dispersion mechanisms in thin graphite-epoxy fiber reinforced composite 
specimens. 

This research is a necessary prelude to further studies on attenuation 
and damage in thin plates or structures of fiber reinforced materials. 





CHAPTER II 


LAMB WAVE EXCITATION 

2.1 TRANSDUCER/SPECIMEN TESTING CONFIGURATIONS 
18 17 21 

Worlton , Victorov , and Firestone have presented methods for 

exciting Lamb waves; the wedge method of inserting Lamb waves of specified 

phase velocities is most common and was adapted for this research. 

Three transducer/specimen configurations were used for testing in this 

research (Figure 22). The configurations where 6*0° and 6=90° are widely 

7 13 18 

used especially in the work of Vary and Williams . Worlton varied 0 

between 0° and 90° as he studied Lamb waves in aluminum and zirconium plates. 

This very last configuration is used to generate and detect Lamb-wave modes 

at a specified angle of incidence. Detecting these modes depends not only 

on the angle of incidence but also on the frequency of the transmitted 

wave, on the material thickness, and on the material properties. 

In all three configurations the wave will be propagating in the 

direction with the parti cal motion in both the x-j and x^ directions. 

Panametric and Dunegan transducers were used in these experiments. 

Their properties are tabulated in Table II on page 6 of Part I. 

Strain gauges were mounted on to the transducer holders to measure 

1 3 

the pressure applied to the transducers. Williams, et.al. defined and 
measured a "saturation pressure" which is required for reproducible coupling 
of transducer to specimen; this pressure was monitored by the strain-gauged 
transducer holders. 
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An Ultragel II ultrasonic couplant was used between the transducer 
and the specimen. A very thin coating of couplant was used where effects 
of thickness of the couplant were neglected. A felt material was used on 
the backside of the transducer holder to insure that energy did not enter 
the transducer holder from the specimen. 

2.2 EXPERIMENTAL SYSTEM 

The experimental system is shown in Figure 23; parts of the Arenberg 
high frequency pulsing system are described below: 

Pulsed Oscillator P.G.650.C 

The P.G.650.C is a high powered pulsed oscillator which provides a 
voltage for the sending transducer and produces a narrow band of frequencies 
which are vital to our measurements since all effects to be measured are 
frequency dependent. The P.G.650.C possesses a wide frequency range from 
.2 MHZ to 190 MHZ with a variable pulse length of 1.5 to 35 usee and 
dele(yed trigger action. It has a refined degree of stability and "fine" 
as well as "coarse" controls on the major variables. It has the flexibility 
of operation as an independent unit or with triggers and modulation provided 
by an external source to synchronize with another system of test equipment. 

Low Frequency Inductance Tuner L.F.T.500 

The L.F.T.500 is designed to be used with the pulsed oscillator P.G.650.C 
using a low frequency current tuned coil. Below 500 Kc, the capacitors 
needed to tune the oscillator would be excessively large; consequently, the 
inductance is varied in each coil by passing a stable D.C. current through 
the output winding on the ferrite cores in a manner which does not interfere 
with the normal operation. Its current range is from 5 Ma to 5A. Its 
calibration linearity is .25%. 


I II II 


I 



FIG. 23 EXPERIMENTAL SYSTEM. 

* The star system was used when frequency 
lower than 500 KC was needed. 

RT Receiving Transducer 
ST Sending Transducer 
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Attenuators ATT 698 

The attenuator furnish the basis for precision measurement of inser- 
tion loss, absorption, etc. in ultrasonic work. The frequency range is 
from O.C. to beyond 100 MHZ. 

2.3 MATERIALS TESTED 

The specimens used were eight-ply graphite-epoxy fiber reinforced 
composites and aluminum (6061,2024) which was used in the last phase of 
this work. Physical properties of the specimens are given in Table X. 

2.4 PROCEDURE 

An ultrasonic wave was transmitted into the specimen by a sending 
transducer at an angle of incidence 0 (see Figure 22). At variable dis- 
tances from the sending transducer in the direction and on the same 
face (except for 0=90°), a receiving transducer was placed to pick up the 
transmitted signal at a receiving angle equal to the angle of incidence. 

The sending transducer input as well as the receiving transducer output 
were displayed on the oscilloscope for real time analysis and frequency 
analysis. Typical input and output signals are shown in Figure 24. 



Specimen 

Number 

Tested 

Ply 

Angle, deg. 

Length, 

L,cm 

Thickness 

h,cm 

Density 

g/cm3 

Ultimate 

Strength 

GPA 

Tensile 

Modulus 

E,GPA 

Composite* 








LVllK 

1 

[0]8 

25.4 

.1384 

1.584 

1.24 

no 

LV33H 

1 

[0/+45/0]S 

25.4 

.1355 

1.355 

.703 

83.4 

LV55 

1 

[+45 ]S 

25.4 

.1340 

1.340 

.181 

20.1 

LV22 

1 

[90]8 

25.3 

.1329 

1.603 

.0672 

8.27 

A1 uminum 


Wi dth 
w,cm 






2024 T4 

2 

1.27/1.57 

25.4 

.159 

2.768 

4.69 

73.1 

2024 T4 

1 

12,7 

25.4 

.254 

2.768 

4.69 

73.1 

2024 T4 

1 

15.88 

25.4 

.318 

2.768 

4.69 

73.1 

6061 TO 

1 

25.4 

8.11 

.081 

2.713 

1.17 

68.95 

6061 T4 

1 

25.4 

16.0 

.159 

2.713 

2.41 

68.95 

6061 T4 

1 

25.4 

11.4 

.229 

2.713 

1.17 

68.95 


TABLE X - PHYSICAL PROPERTIES OF THE SPECIMENS 

* after Vary® 
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FIG. 24 


TYPICAL INPUT AND OUTPUT SIGNALS AND THEIR 
CORRESPONDING FREQUENCIES FOR LVIIK. 
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CHAPTER III 
EXPERIMENTAL RESULTS 


Experimental results are presented in this chapter for the 
dispersion phenomenon for the four composite specimens LVllK, LV33H, 

LV55, and LV22 in the case of normal incidence i.e. 0=0° and in the 
case of axial incidence i.e. 0=90°. A discussion and list of advantages 
of using plastic v/edges in the case where 0°<9<90° will be presented. 

And finally the results of the testing of the Arenberg ultrasonic system 
in detecting Lamb-wave in aluminum (2024,6061) in the axial incidence case 
will be presented. 

3.1 DISPERSION IN THE CASE OF NORMAL INCIDENCE 

Wave speed measurements were recorded as waves of different frequencies 
(f = .15 - 3.0 MHZ) were traversing the four composite specimens. A graph 
ms plotted showing speed versus frequency for the four composite specimens 
under studyin Figure 25. The graph shows that the wave velocity is frequency 
dependent for all four specimens. 

The bulk velocities for LVllK and LV22 are approximately 850 cm/ms and 
100 cm/ms respectively. By comparing these values to the velocities recorded 
and presented in Figure 25 for LVllK and LV22, one can see that the corresponding 
wave velocities were higher than the bulk velocity for the same material. 

This is explained by the prese,nce of Lamb waves. For the specimens LVllK and 
LV22 the dispersion data revealed the existence of higher modes of wave propa- 
gation. In order to define these higher modes plexiglass wedges (O°<0<9O°) 
v;ere used. 
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Spend C, cm/ms 


1100 

1070 

1010 

1000 

950 

910 

880 

850 


AA 


A LV11K 
□ LV33H 
■ LV55 
• LV22 


IK 


750 

700 


hn 


□ 


□ 


□ 




550 


□ 


400 

330 

300 

270 

200 

170 


100 


“ (C^)lV22 

Frequency f, MHZ 


.15 


Fig. 25 Wave speed vs. frequency for composites with 
normal incidence transducer arrangement. 
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3-2 THE USE OF PLASTIC WEDGES FOR 0® < 0 <90® 

22 

Plona, Behravesh, and Mayer have shown that the phase velocity 
of Lamb waves may be controlled by setting the angle of incidence of 
an immersing longitudinal wave as shown in Figure 22a; i.e. by vaiying 6. 
This may be done by immersing the transducer in liquid, in air, or by 
using a plastic wedge. 10, 15 and 30 degrees plastic wedges were made 
and tried in this research, but evidence is presented only for the 15® 
wedge. Great difficulty was encountered in attempting to hold the wedges 
properly. Wave velocity data was not collected, but qualitative evidence 
of the value of using wedges will be noted. 

Figure 26 compares three oscilloscope time and frequency traces; 
the material is LVllK and the frequency is approximately 1 MHZ. The 
input signal is always the wave in the upper left hand comer of the trace. 
The top trace shows the output for normal incidence (0=0®) where the effects 
of constructive and destructive interference are seen due to the multiple 
reflections across the thickness of the specimens. The output trace appears 
as a modulated wave and is clearly quite different in form or shape from 
the input wave. It is also probable that multiple modes and end reflections 
are present. The middle trace shows the output for axial incidence (0=90®) 
where the effects of multiple modes and end reflections are present. The 
output is once again quite different in form and shape from the input. The 
bottom trace shows the output when a 15® plastic wedge is used. In this 
case the output wave has almost exactly the same shape and form as the input 
wave, i.e., the time duration is 30 usee and the amplitude is constant. It 
is believed that the wedge method eliminates the constructive/destructive 
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FIG. 26 - TRACES SHOWING WEDGE EFFECT VERSUS KONWEDGE EFFECT 
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interferences effect, the multiple mode problem, and the end reflection 
waves. 

Figure 27 shows how the wedge eliminated a problem which frequently 
occurred for frequencies greater than 1 MHZ; electromagnetic waves were 
often transmitted when the transducer was in direct contact with the 
specimen; the material is LVllK and the frequency is approximately 2.5 MHZ. 
The top trace shows the output for axial incidence (0=90°) where two 
different waves are present, one wave arriving instantenously and the other 
wave arriving later. The first arrival is an electromagnetic wave whose 
speed is obviously too great to measure. The second arrival is the mechani- 
cal wave response for the input wave. The output signal shown in the lower 
trace, where a 15° wedge was used, shows no evidence of an electromagnetic 
wave. 


3-3 AXIAL INCIDENCE WAVES, 0=90° 

This section describes how the lowest order symmetric Lamb waves are 

propagated along the axial direction (x^) of the four composite specimens. 

In 'Figure 28 C versus fh was plotted for all four specimens, where C is the 

velocity of the stress wave, f is the frequency of the propagating wave, 

and h is the thickness of the specimen. 

Figure 28 only presents a small portion of the total possible dispersion 

22 

curves for this material. Habeger, Mann, and Baum have developed a theory 
which may be used to construct dispersion curves for the 0° degree (LVllK) 
and the 90° (LV22) unidirectional laminates. The authors of this report are 
not able to construct these curves at this time because we do not have the 
specimens necessary to conduct tests to measure material constants needed for 



64 




b. OBLIQUE INCIDENCE 
1 5 WEDGE 


FIG. 27 - TRACES SHOWING THE WEDGE EFFECT IN ELIMINATING 
ELECTROMAGNETIC WAVES . 


fh. MHZ-cm 


Fig. 28 C vs. fh for axial incidence wave, 0 

A LV11K 
□ LV33H 
■ LV55 
• LV22 
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the theory. It is our belief that the data presented in Figure 28 would 

lie on the initial portion of the lowest order symmetric mode. Note the 

largest value of fh plotted on the abscissa is 0.4 MHZ-cm. The thickness 

for these specimens is approximately 0.13 cm and the maximum available 

frequency for testing is 3.0 MHZ; the frequency is limited by the available 

transducers. The data for LVllK shows a slight decrease as the product fh 

increases; this is consistent with the results reported by Habeger, Mann, 

22 

and Baum . A limited amount of data is available for LV33H, LV55, and 
LV22 because of great energy dissipation with increasing frequencies over 
the long length of the specimen. 

3-4 LAMB WAVE PROPAGATION IN ALUMINUM 

Lamb wave propagation in aluminum was studied to demonstrate that our 

experimental system and methodology was working and was in agreement with 

the analytical results^^, in particular we studied axial incidence (0=90°) 

waves and compared them with theory. 

The data was plotted on C versus fh curves from the work done by 
1 8 

Worlton for aluminum plate, and is presented in Figures 29 and 30. 

Aluminum specimens of five different thicknesses were used in order to 
generate more data points. 

It is not possible to draw physically correct curves through the data 
points without the aid of an analytical theory. If one removes the families 
of curves (modes), the experimental data appears randomly scattered. 

Because the experimental results reported for aluminum using the 
Arenberg systems are in agreement with existing theory and experiment, it 
may be concluded that the system is functioning correctly and that Lamb wave 
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fhxlO® cy-in/sec-^ 


Fig. 29 Lamb wave propagation in aluminum 
(6061). 

• h = .09 in 
o h = .032 in 
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propagation exists in specimens where thickness is of the order of the 
composite specimen thickness and hence it is correct to assume Lamb waves 
are propagating in the composites. It is also apparent that an analytical 
theory must be used to develop dispersion curves for composites before the 
composite data may be plotted. 
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CHAPTER IV 
CONCLUSION 


The data collected for the aluminum and composite specimens 
clearly demonstrate that Lamb wave propagation is the mode of 
propagation in existence for all angles of incidence, 

A theory is available which may be applied to homogeneous 
orthotropic materials which, for this study, were the 0° degree and 
the 90° degree graphite-epoxy fiber reinforced composite. Additional 
specimens and testing are needed to apply the theory. The [0/-45/0]s 
and the [-45 ]s graphite-epoxy fiber reinforced composite which are 
not homogeneous orthotropic materials will require a more advanced 
theory. 

Plastic wedges appear to be a feasible way of sending and 
receiving Lamb waves at oblique incidence, but more work must be 
done on constructing a device to hold transducers and wedges at 
arbitrary angles. 
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